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Rate  constants  and  product  ion  branching  fractions  for  the  gas-phase  reactions  of  02  and  O3  with  the  anions 
(a)  P02C1-,  (b)  POCl3-,  (c)  POCl2-,  and  (d)  P02C12-  were  measured  in  a  selected-ion  flow  tube  (SIFT).  The 
kinetics  were  measured  at  temperatures  of  163—400  K  and  a  He  pressure  of  0.4  Torr.  Only  P02C1~  reacts 
with  02  to  a  measurable  extent,  having  £(163—400  K)  =  1.1  x  10-8(77K)~10  cm3  molecule-1  s-1,  while  03 
reacts  with  all  of  the  anions  except  P02C12-.  The  fitted  rate  constant  expressions  for  the  03  reaction  with 
anions  a-c  are  as  follows:  £a(163-400  K)  =  3.5  x  10-6(77K)-L6,  £b(163— 400  K)  =  4.0  x  10-7(27K)-1-2, 
and  £c(163— 400  K)  =  3.7  x  10-7(77K)-L4  cm3  molecule-1  s-1.  Calculations  were  performed  at  the  G3  level 
of  theory  to  obtain  optimized  geometries,  energies,  and  electron  affinities  (EAs)  of  the  reactant  and  product 
species,  as  well  as  to  determine  the  reaction  thermochemistry  to  help  understand  the  experimental  results. 
The  PO.vCly-  anions  that  have  lower  electron  binding  energies  (eBE)  and  higher  spin  multiplicities  are  more 
reactive.  The  doublets  are  more  labile  than  the  singlets.  How  the  extra  electron  density  is  distributed  in  the 
anion  does  not  predict  the  observed  reactivity  of  the  ion.  The  reactions  of  P02C1-  with  02  and  O3  yield 
predominantly  P03-  and  PO4-.  The  reaction  of  POCl2-  with  03  yields  mostly  Cl-  and  P02C12-,  while  the 
POCI3-  reaction  with  03  yields  mostly  03-  and  P02C12-. 


Introduction 

The  toxicity,  reactivity,  and  decomposition  pathways  of 
phosphorus-containing  compounds  are  particularly  interesting 
for  their  role  in  the  use,  production,  incineration,  and  accidental 
emissions  of  chemical  warfare  nerve  agents,  pesticides,  plasti¬ 
cizers,  flame  retardants,  and  hydraulic  fluids.1-6  Concerns  arise 
in  the  contamination  of  plasma  by  these  compounds  in  coal- 
fired  magnetohydrodynamic  electric  generators.7  Also,  flame 
tests  intimate  that  the  addition  of  a  few  percent  of  either 
phosphine  (PH3)  or  oxyphosphorus  compounds  to  hydrogen- 
fueled  scramjet  combustors  yields  a  measurable  increase  of  the 
flameholding  stability  and  output  thrust  in  the  supersonic 
expansion  because  of  enhanced  H,  O,  and  OH  recombination 
kinetics.8-12  Furthermore,  evidence  suggests  that  combustion 
is  sensitive  to  and  can  be  enhanced  by  ion— molecule  chemis¬ 
try.13  Consequently,  the  reactions  of  certain  neutral  and  anionic 
phosphorus-containing  compounds  with  common  gases  have 
been  studied,14-18  in  part  to  elucidate  pathways  for  P03- 
production.18’19 

To  address  these  issues,  kinetics  measurements  for  the 
reactions  of  P02C1-,  P02C12-,  POCl2-,  and  POCl3-  anions  with 
the  oxidants  02  and  O3  (ozone)  have  been  made.  Ozone  is 
interesting  because  it  can  attach  electrons  and  it  is  a  very  strong 
oxidant,  while  reactions  with  group  V  compounds  producing 
species  such  as  P02  can  be  chemiluminescent.20-22  Rate 
constants  and  product  branching  fractions  have  been  measured 
from  1 63  to  400  K  for  the  POjCly  ions  mentioned  above  and 
these  data  are  compared  with  previous  work  on  these  systems. 
However,  a  dearth  of  studies  of  the  ion  chemistry  of  these 
systems  exists.16’23-25  The  results  of  calculations  at  the  G3  level 
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of  theory  to  obtain  optimized  geometries,  0  K  total  energies, 
298  K  enthalpies,  and  electron  affinities  (EA)  are  also  included 
for  the  relevant  reactant  and  product  molecules.  These  results 
have  been  utilized  to  determine  the  reaction  thermochemistry 
and  to  help  understand  the  experimental  data. 

Experimental  Section 

The  AFRL  selected  ion  flow  tube  (SIFT)26’27  was  used  to 
make  the  measurements.  The  instrument  has  been  discussed  in 
detail  elsewhere  and  is  only  briefly  described  as  it  pertains  to 
the  current  experiments.  The  ions  were  formed  in  a  remote 
moderate  pressure  electron  impact  ion  source  at  about  0.1  to  1 
Torr.  P02C1-,  P02C12-,  POCl2-,  and  POCl3-  were  formed  from 
POCI3  (Aldrich,  99%).  The  POCl3  was  used  as  obtained  from 
the  manufacturer  except  for  perfonning  several  freeze— pump- 
thaw  steps  to  remove  condensables.  Reactant  ions  were  selected 
with  a  quadrupole  mass  filter  and  then  injected  into  a  fast  flow 
of  helium  buffer  gas  (AGA,  99.997%)  in  a  1-m  long  flow  tube 
using  a  Venturi  inlet,  giving  a  flow  tube  pressure  of  0.4  Torr. 
The  He  buffer  gas  was  passed  through  a  liquid  nitrogen  cooled 
sieve  trap  in  order  to  remove  water  vapor  before  entering  the 
flow  tube.  The  reactant  ions  were  then  thermally  equilibrated 
and  proceeded  downstream  past  an  inlet  used  to  introduce  the 
neutral  reactant.  Both  02  (Massachusetts  Oxygen,  99.999%)  and 
03  were  used  alternately  in  this  study.  The  ozone  was  produced 
by  interfacing  an  Orec  O3V-0  ozonator  to  the  SIFT  apparatus 
as  described  in  detail  in  a  previous  work.14  This  technique 
yielded  a  reactant  flow  of  ca.  5%  03  in  02.  At  the  end  of  the 
reaction  tube,  a  portion  of  the  gas  was  sampled  through  an 
aperture  in  a  blunt  nose  cone  into  a  second  quadrupole,  and  the 
mass-filtered  ions  were  counted  using  a  particle  multiplier. 

The  rate  constants  were  obtained  from  the  pseudo-first-order 
attenuation  of  the  reactant  ion  signal  that  depends  on  the  reaction 
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TABLE  1:  G3°  Total  Energies,  Enthalpies,  and  Electron 
Affinities  (EA),  All  Units  except  EA  in  Hartrees 


species 

total 

energy 

(OK) 

enthalpy 
(298  K) 

EA5  (eV) 

calcd 

EA® 

(eV) 

Ph  cl, 

0  b 

-75.03099 

-75.08014 

-75.02863 

-75.07778 

1.338  (1.461)'' 

o2(3z-g) 

o2-(2ngi) 

-150.24821 

-150.26294 

-150.24490 

-150.25963 

0.401  (0.451)® 

03(C2„,>A,) 

o3-(c2„,2b,) 

-225.31826 

-225.39754 

-225.31437 

-225.39357 

2.157  (2.1028/ 

C1(2P) 

cr(‘S) 

-459.99096 

-460.12360 

-459.98860 

-460.12124 

3.609  (3.6127/ 

chCz^) 

cir(2z\) 

-920.07127 

-920.16188 

-920.06777 

-920.15805 

2.465  (2.38/ 

0  0 
°  0 

-535.12026 

-535.20556 

-535.11685 

-535.20200 

2.321  (2.276/ 

C102(C2„,2B,) 

ao2-(C21„1A1) 

-610.24282 

-610.32562 

-610.23871 

-610.32123 

2.245  (2.140/ 

C120(CV>Ai) 

C120"(G,2A7' 

-995.16647 

-995.25362 

-995.16208 

-995.24808 

2.37  (>2.2/ 

P02(C2w2A,) 

P02-(C2i,,'A') 

-491.59301 

-491.72030 

-491.58890 

-491.71624 

3.464  (3.42/ 
3.7/  3.6/ 

P03(C2,.,2B2) 

PO  3-(D3;„1A1) 

-566.77102 

-566.96044 

-566.77622 

-566.95590 

5.155(4.95/ 

P04(C2„,2A2) 

po4-(c2i„>a,) 

-641.884035 

-642.070234 

-641.878008 

-642.064763 

5.067" 

POC12(Cs,2A') 

POC12-(Cs,1A') 

-1336.5486 

-1336.6857 

-1336.5429 

-1336.6795 

3.732  (3.8/ 

3.71 

POCl3(CWA,) 

POC13-(Cs,2A') 

-1796.67736 

-1796.73562 

-1796.67057 

-1796.72744 

1.586(1.4/ 

1.50 

POjCKGJA,) 

P02C1-(C.„2A') 

-951.710965 

-951.789794 

-951.705948 

-951.784401 

2.145 

P02C12(C,,2A) 

po2ci2-(C2„,‘a,) 

-1411.72046 

-1411.93316 

-1411.71449 

-1411.92673 

5.788 

P03C1(C(,'A) 

P03C1-(C„2A) 

-1026.81474 

-1026.97642 

-1026.80895 

-1026.97046 

4.399 

“  The  Gaussian  2003  program  package30  was  used  to  perform  the 
calculations.  b  The  numbers  in  parentheses  are  experimental  results  from 
previous  work.  c  Calculated  EA  using  the  G2  method.47  d  Neumark 
and  Lykke.48  e  Travers  et  al.49  f  Drzaic  et  al.50  *  Berzinsh  et  al.SI 
h  Gilles  et  al.52  1  Wecker  et  al.53  >  The  structure  of  this  ion  (see  Table 
2)  has  a  long  Cl— O  bond  that  indicates  a  more  electrostatic  interaction 
which  may  introduce  inaccuracies  in  the  Hartree— Fock  zero  point 
energy  (ZPE)  value  since  electron  correlation  is  not  taken  into  account. 
To  correct  for  this,  an  MP2  calculation  was  performed  and  the  ZPE 
was  found  to  be  28  meV  greater.  Thus,  a  more  accurate  estimate  of 
the  EA  would  be  lower  by  this  value.  k  Xu  et  al.54  1  Miller.55  m  Lohr.22 
"  Wang  and  Wang.19  0  EAfPO/  was  calculated  using  the  G3B3  method 
(see  text).  p  Mathur  et  al.23 

time  and  the  concentration  of  the  neutral  reactant,  the  latter  of 
which  is  varied  during  a  measurement.  The  rate  constants  have 
relative  uncertainties  of  ±15%  and  absolute  uncertainties  of 
±25%.26  For  the  product  distributions,  the  SIFT  was  operated 
under  conditions  of  modest  attenuation  of  the  reactant  ion  signal 
and  the  product  branching  fractions  were  extrapolated  to  a 
condition  of  zero  neutral  reactant  concentration.  Thus,  the  effects 
of  secondary  chemistry  are  accounted  for.  The  mass  resolution 
of  the  downstream  mass  filter  has  also  been  adjusted  to  minimize 
mass  discrimination  during  the  branching  fraction  measurements. 
The  branching  fractions  have  relative  uncertainties  of  ±10% 
of  the  major  product  peaks.28 

Because  of  the  unavoidable  presence  of  O2  during  the  O3 
experiments,  corrections  to  the  branching  fractions  for  the 
PO2CI  reaction  with  O3  were  made  by  subtracting  the  product 
counts  from  the  02  reaction  with  P02C1~.  When  the  POCl3_  ± 
03  system  was  studied,  the  breakup  of  roughly  5%  of  POCl3_ 


TABLE  2:  Bond  Lengths  (in  A)  and  Angles  (in  deg)  for 
Relevant  Species  Obtained  at  the  G3  Level  of  Theory, 
Where  the  Superscripts  Distinguish  the  O  and  Cl  Atoms 


parameter 

po2cr 

PO2CI2- 

PO2CI2  P02C1 

par  P03C1 

por  ci2o- 

r(P-O) 

1.509 

1.491 

1.630 

1.477 

1.666°  1.470* 

1.508 

r(P— O2) 

1.477 

1.666  1.620 

r(P— O3) 

1.496  1.620 

r(P-Cl') 

2.270 

2.115 

2.010 

1.998 

1.987 

tfCP-O) 

1.700 

r(Cl2— O) 

2.134 

Z(0‘-P-C1‘) 

106.28 

106.8 

98.8 

112.7 

114.8 

Z(02-P-C1>) 

108.2 

ZfO'-P-O2) 

125.95 

127.5 

119.8 

134.7® 

57.7°  126.7 

120.0 

Z(02-P-03) 

114.6  61.2 

120.0 

Z(03-P-04) 

123.2 

Z(C1‘-P-C12) 

98.3 

106.9 

Z(CP— O— Cl2) 

116.3 

dihedfC/O2/ 

136.3 

102.0 

180.0 

147.6 

dihed(0',Cl2/ 

110.5 

dihedfO'.Cl1/ 

127.8 

125.7 

142.7 

dihed(0',03/  104.6 

dihed(0',02)*  64.1 

“  PO4-  has  two  double  bonds  and  two  single  bonds.  O  atoms  labeled 
1  and  2  are  single-bonded  while  3  and  4  are  double-bonded  to  P. b  O 
atom  1  is  double  bonded  to  P  while  the  other  two  share  a  bond.  c  The 
O— P— O  angle  agrees  with  the  bounded  estimate,  130 — 140°.4S  rfP— 
Cl1  is  the  reference  axis. e  P-O1  is  the  reference  axis.  7P— O2  is  the 
reference  axis.  g  P—O3  is  the  reference  axis. 

to  POCI2  upon  injection  into  the  flow  tube  was  corrected  by 
a  similar  subtraction  procedure  because  POCl2~  is  also  reactive 
with  03.  Another  breakup  product  is  Cl~;  however,  a  back¬ 
ground  subtraction  could  also  be  used  as  Cl"  is  unreactive  with 
02  and  03.29  These  observations  are  consistent  with  previous 
studies  of  both  dissociative  and  nondissociative  electron  attach¬ 
ment  to  POCI3,  where  POCl2“  and  Cl~  are  the  major  and  minor 
dissociative  products,  respectively.15,16’25  In  all  cases,  the  product 
branching  ratios  were  renormalized  after  making  background 
corrections. 

Computational  Method 

The  Gaussian  03W  program  package30  was  used  to  perform 
calculations  on  the  reactant  and  product  species  involved  using 
the  G3  theoretical  method.31  The  wave  function  stability  of  the 
minimum  energy  structures  from  the  Hartree— Fock  (HF) 
optimization  step  of  G3  at  the  6-31G(d)  level  have  been  verified 
for  the  neutral  and  ionic  species,  with  the  exception  of  PO4 
neutral  as  discussed  later.  The  G3  0  K  energies,  298  K 
enthalpies,  and  EAs  are  compared  with  the  literature  values 
whenever  possible.  These  data  are  given  in  Table  1  for  the 
optimized  structures  that  have  the  parameters  listed  in  Table  2. 
The  G3  enthalpies  at  298  K  have  also  been  used  to  calculate 
the  heat  of  reaction  at  298  K,  A//rxn,  for  all  of  the  pathways 
discussed  in  the  following  sections.  The  average  absolute 
deviation  in  the  G3  results  is  ±0.041  eV  (0.94  kcal  mol-1).31 
Additionally,  a  natural  population  analysis32  (NPA)  was  used 
to  ascertain  the  location  of  the  “extra”  electron  density  in  the 
anions;  geometries  were  optimized  using  MP2(Full)  at  the 
6-3  ll±G(d)  level. 

Experimental  Results 

Rate  constants  and  branching  fractions  for  the  reactions  of 
P02C1“,  P02C12±  POCl2_,  and  POCl3_  with  02  and  O3  were 
measured  from  163  to  400  K  at  a  He  pressure  of  0.4  Torr.  Only 
PO7CI  reacts  with  02  to  a  measurable  extent,  while  P02C12“ 
was  found  to  be  unreactive  with  both  02  and  03.  This 
observation  is  in  agreement  with  earlier  work  where  P02C1?_, 
POCl2~,  and  POCl3_  were  found  to  be  unreactive  with  02. 17,25 
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Figure  1.  Rate  constants  in  cm3  s-1  for  the  reactions  of  P02C1-, 
POCl2-,  and  POCE-  with  O2  and  O3  as  a  function  of  temperature  in  K 
as  measured  in  a  selected  ion  flow  tube  (SIFT). 

Moreover,  this  lack  of  reactivity  is  not  surprising  as  the  EA  of 
O2  is  smaller  than  that  of  POCI3  (see  Table  1).  The  SIFT 
apparatus  is  capable  of  accurate  measurement  of  rate  constants 
approximately  as  low  as  5  x  10-12  cm3  s-1.  Consequently,  the 
rate  constants  for  the  reactions  of  PO2CI2-,  POCE-,  and  POCI3- 
with  O2  have  a  value  <5  x  10~12  cm3  s-1.  A  similar  upper 
limit  can  be  placed  on  the  rate  constant  for  the  reaction  of 
P02C12-  with  O3  as  it  is  unreactive  as  well,  since  no  products 
are  formed. 

The  rate  constants  for  the  reaction  of  POiQ  with  02  given 
by  eq  1 

2P02CF  +  302  —  *P03“  +  2C10  A//rxn  =  —1.18  eV  (la) 
'P04~  +  2C1  A Han  =  -0.48  eV  (lb) 

are  shown  in  Figure  1.  The  superscripts  in  eq  1  and  all 
subsequent  reaction  equations  denote  the  spin  multiplicity  and 
only  the  spin-allowed,  exothermic  channels  are  given.  The 
energetics  of  the  reactants  and  products  calculated  using  G3 
theory  for  all  the  listed  reactions  are  also  shown  in  Table  1. 
For  puiposes  of  efficiency,  the  rate  constant  data  for  all  of  the 
reactions  studied  are  shown  in  Figure  1.  A  best  fit  to  the  rate 
constants  for  reaction  1  in  Figure  1  yields  the  expression  A:(  1 63 — 
400  K)  =  1.1  x  10-8(77K)-10  cm3  molecule-1  s-1.  This  case 
is  the  only  one  where  oxygen  reacts  to  a  measurable  extent  with 
PCE-Clv-  in  the  SIFT.  Clearly,  the  rate  constants  decrease  with 
increasing  temperature,  indicative  of  a  reaction  involving  either 
an  intermediate  collision  complex  or  transfer  of  an  electron. 

The  product  ion  branching  fractions  for  the  products  P03- 
and  PO4-  are  listed  in  Table  3.  The  fraction  of  PO4-  is 
somewhat  greater  than  that  of  PO3-  at  all  temperatures  and 
increases  with  increasing  temperature. 

The  rate  constants  for  the  reaction  of  P02C1-  with  O3  given 
in  eq  2 

2P02C1_  +  '03  —  'P03“  +  2C102  M/rxn  =  -2.61  eV  (2a) 


TABLE  3:  Rate  Constants"  and  Ion  Product  Branching 
Fractions  for  the  Reactions  of  PO2CI-,  PO2CI2-,  POCE-, 


and  POCI3- 

with  02  and  O3* 

reaction 

163  K 

294  K 

400  K 

PO2CI-  +  02 

0.60  [6.0] 

0.33  [6.0] 

0.24[6.0] 

— - 

P03-  +  CIO 

0.49 

0.42 

0.41 

PO4-  +  Cl 

0.51 

0.58 

0.59 

P02C1-  +  03 

1 1  [9.9] 

5.0  [8.9] 

2.6  [8.6] 

P03-  +  C102 

P03-  +  Cl  +  02 

0.62 

0.65 

0.66 

PO4-  +  CIO 

0.34 

0.34 

0.33 

Cl"  +  po3  +  o2 

0.04 

0.01 

0.01 

pocir  +  03 

3.3  [9.7] 

1.1  [8.6] 

1.0  [8.3] 

* 

Cl-  +  PO2CI  +  02 

Cl-  +  P03  +  CIO 

0.71 

0.79 

0.60 

PO2CI2-  +  02 

0.26 

0.20 

0.35 

PO3-  +  CIO  +  Cl 

P03-  +  cl2o 

0.02 

0.01 

0.03 

P04-  +  Cl2 

0.01 

<0.005 

0.02 

POCI3-  +  O3 

10  [9.3] 

4.0  [8.3] 

3.6  [8.1] 

-*■ 

03-  +  POCE 

0.50 

0.55 

0.55 

P02C12-  +  C102 
P02C12-  +  Cl  +  02 

0.44 

0.42 

0.44 

Cl-  +  PO2CI2  +  o2 

Cl-  +  P02  +  Cl2  +  02 
Cl-  +  PO3CI  +  CIO 

0.06 

0.03 

0.01 

"Rate  constants  in  units  of  10-10  cm3  s-1;  bracketed  numbers 
represent  calculated  collisional  rate  constants  obtained  by  the  param¬ 
etrized  trajectory  approach  of  Su  and  Chesnavich.56'S7  b  Measurements 
of  the  rate  constants  of  P02CE-,  POCl2~,  and  POCI3-  with  02  and  of 
PO2CI2-  with  O3  yielded  upper  limits  of  5  x  10-12  cm3  s-1. 


are  also  shown  in  Figure  1,  while  the  product  branching  fractions 
are  listed  in  Table  3.  A  best  fit  to  the  rate  constants  gives  the 
expression  £(163—400  K)  =  3.5  x  10-6(77K)-1-6  cm3  molecule-1 
s-1  for  reaction  2.  Clearly,  O3  reacts  about  an  order  of  magnitude 
faster  than  02.  The  same  decrease  with  increasing  temperature 
is  also  seen  with  ozone.  Unlike  the  previous  reaction,  the  major 
product  ion  is  PO3-  and  its  fraction  is  about  twice  as  great  than 
that  of  PO4-.  Also,  a  small  amount  of  Cl-  is  produced;  however, 
it  is  less  than  5%  of  the  total  products  and  it  also  decreases 
with  increasing  temperature. 

The  rate  constants  for  the  reaction  of  POCl2~  with  03  given 
by  eq  3 


‘POCl2”  +  '03  —  'cr  +  'P02C1  +  '0,(a'Ag) 

Aft,™  —  —1.15  eV  (3a) 


'cr  +  2po3  +  2cio 

'P02C12-  +  '02(a'Ag) 
'P03_  +  2C10  +  2C1 
'P03“  +  'ci2o 
'p<V  +  'ci2 


AF/^n  —  —0.557  eV  (3b) 
AFfjxn  —  —3.86  eV33  (3c) 
A//™  =-1.84  eV  (3d) 
Atfrxn  =  -3.34eV  (3e) 

A =  -3.60  eV  (3f) 


'P03_  +  2Cl  +  302 
'P03~  +  2C1  +  ‘02(a'Ag) 
‘P04~  +  2C10 

'cr  +  2po3  +  3o2 
'cr  +  2P03  +  '02(a'Ag) 


A//rxn  =  -2.47eV  (2b) 
A7/rxn  =  -1.49  eV  (2c) 
A7frxn  =  — 2.08  eV  (2d) 
A//rxn=-1.19eV  (2e) 
A/4n  =  -0.21eV  (2f) 


are  well  fitted  by  the  expression  £(163—400  K)  = 
3.7  x  10-7(77K)-1-4  cm3  molecule-1  s-1  and  are  somewhat  less 
than  the  rate  constants  for  the  P02C1-  reaction.  Nevertheless, 
the  rate  constants  demonstrate  the  same  negative  trend  with 
increasing  temperature.  The  product  ion  branching  fractions  in 
Table  3  show  that  Cl-  is  by  far  the  major  product,  while  the 
P02C12-  yield  is  roughly  25%  and  increases  to  35%  at  400  K 
at  the  expense  of  Cl-.  Trace  amounts  of  PO3-  and  PO4-  were 
also  observed. 
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TABLE  4:  Natural  Population  Analysis  of  the  PO*Cl,, 
Species" 


charge  associated  with  each  atom 

P 

O1 

O2 

Cl1 

Cl2 

Cl3 

pocb 

1.547 

-0.95 

-0.199 

-0.199 

-0.199 

pocir 

1.399 

-1.012 

-0.404 

-0.492 

-0.492 

extra  e~ 

0.15 

0.06 

0.21 

0.29 

0.29 

poch 

1.405 

-0.929 

-0.238 

-0.238 

poet" 

1.117 

-1.022 

-0.547 

-0.547 

extra  e~ 

0.29 

0.09 

0.31 

0.31 

P0,C124 

1.777 

-0.435 

-0.944 

-0.199 

-0.199 

P02C12" 

1.952 

-1.063 

-1.063 

-0.413 

-0.413 

extra  e~ 

-0.175 

0.628 

0.119 

0.214 

0.214 

PO2CI 

1.980 

-0.883 

-0.883 

-0.214 

poor 

1.630 

-1.048 

-1.048 

-0.533 

extra  e~ 

0.350 

0.165 

0.165 

0.319 

a  Superscripts  used  to  label  different  O  and  Cl  atoms.  b  Here,  the 
charges  on  the  two  O  atoms  are  not  equivalent  since  the  P—O  bond 
lengths  are  not  equal;  see  Table  2. 

Finally,  the  reaction  of  POCI3-  with  O3  given  by  eq  4 

2P0C13~  +  ‘03  — -  203_  +  'POCl3  M/ran  =  -0.61  eV  (4a) 

‘P02C12~  +  2C102  A/7rxn  =  -3.36  eV  (4b) 

'P02C12~  +  2C1  + 30,  A//rxn  = -3.22  eV  (4c) 

’PO.Cl,"  +  2C1  +  'o2(a'Ag)  A Hnn  =  -2.24  eV  (4d) 

'cr  +  2P02C12  +  302  A//rxn  =  —  1 .06  eV  (4e) 

'CP  +  2P02C12  +  ‘02(a'Ag)  A Hnn  =  -0.08  eV  (4f) 

'cr  + 'P03C1  +  2C10  A//rxn  =  — 0.14  eV  (4g) 

is  as  rapid  as  that  of  P02C1"  with  O3.  The  best  fit  expression  to 
the  rate  constants  is  &(  163 —400  K)  =  4.0  x  10"7(77K)"1,2  cm3 
molecule"1  s"1  for  reaction  4.  O3"  is  the  major  product,  resulting 
from  direct  electron  transfer,  while  the  P02C12  product  yield 
is  somewhat  less.  Another  similarity  with  the  P02C1"  reaction 
is  the  presence  of  trace  amounts  of  Cl"  product  ions  which  also 
decrease  with  increasing  temperature. 

Computational  Results 

Table  1  shows  the  G3  calculated  0  K  energies,  298  K 
enthalpies,  and  EAs  of  the  relevant  species.  The  calculated  EAs 
and  the  available  experimental  values  from  earlier  work  are  in 
excellent  agreement.  The  G3(B3LYP)  method  was  used  to 
calculate  EA(PC>4)  because  of  a  problem  in  maintaining  the 
lowest-energy  wave  function  throughout  the  HF  and  Moller— 
Plesset  (MP)  steps  of  an  attempted  G3  calculation  for  neutral 
PO4.  However,  the  G3(B3LYP)  method  is  of  comparable 
accuracy  to  the  G3  method.34 

Table  2  contains  the  optimized  geometries  of  some  of  the 
reactant  and  product  anions.  Examining  the  data  in  Table  2 
shows  that  the  POxCl>,"  anions  studied  in  the  SIFT  all  have 
nonplanar  minimum  energy  structures.  PO2CI"  and  POCl2”  both 
have  umbrella  shaped  geometries,  while  POCI3"  and  PO3CI  are 
tetrahedral.  Adding  the  electron  to  the  PO2CI  neutral  shown  in 
Table  2  causes  the  trigonal  planar  structure  to  distort  to  the 
umbrella  shape.  The  P02C12"  ion  has  the  two  O  atoms  on  the 
same  side  of  the  P  atom  arranged  trans  to  the  two  Cl  atoms  on 
the  opposite  side. 

The  results  of  the  NPA  computations  given  in  Table  4  show 
that  the  electron  densities  of  both  the  neutral  and  charged  species 
are  concentrated  around  the  more  electronegative  atoms,  Cl  and 


O.  Interestingly,  most  of  the  extra  electron  density  from  adding 
an  electron  to  form  the  anion  resides  on  the  Cl  atom(s),  whereas 
the  O  atom(s)  receive  a  much  smaller  portion.  This  was  also 
observed  in  the  NPA  calculation  from  the  electron  attachment 
study  involving  POCI3  and  PCI3.16  The  behavior  of  P02C12"  is 
different  than  that  of  the  other  anions  in  that  the  O  atoms  are 
not  equivalent  and  do  get  a  considerable  portion  of  the  extra 
density.  Moreover,  the  P  atom  of  PO2CI2  loses  electron  density 
upon  attachment. 

Discussion 

It  is  clear  that  ozone  is  more  reactive  than  oxygen,  despite 
the  fact  that  ground-state  ozone  has  a  singlet  multiplicity  while 
ground-state  oxygen  is  a  triplet.  These  results  are  consistent  with 
the  strong  oxidizing  ability  of  ozone.35  The  negative  temperature 
of  the  rate  constants  indicates  that  either  a  collision  complex 
or  a  transfer  of  electrons  is  involved  in  the  reactions.  Ozone 
reactions  with  olefins  can  proceed  via  channels  utilizing  five- 
membered-ring  ozonide  intermediates.35,36  If  a  similar  mecha¬ 
nism  is  involved  in  these  reactions,  the  ozone  could  add  to  the 
double  bond  connecting  P  and  O. 

Alternatively,  a  four-membered  collision  complex  may  be 
formed  by  an  O3  attack  on  the  electropositive  P  atom  of  the 
P0XC1V";37  the  approach  would  be  end-on  since  the  terminal  O 
atoms  of  O3  are  more  negatively  charged  than  the  middle  atom. 
The  NPA  results  in  Table  4  confirm  that  the  P  atom  is  more 
positively  charged.  The  more  positive  middle  O  atom  of  O3 
would  be  attracted  to  the  more  negatively  charged  Cl  atoms  in 
the  collision  complex.  The  spin-allowed,  exothermic  product 
channels  are  listed  in  Table  3  and  reactions  1-4  may  well 
produce  ClmO„  neutral  fragments  or  O2  (a  ’A),  indicating  that 
rearrangements  could  subsequently  occur  during  the  reaction 
that  would  result  in  the  observed  product  ions.  A  similar 
mechanism  occurs  in  the  reaction  of  alkyl  cations  with  O3.29  In 
those  systems,  the  intermediate  complex  forms  with  the  tenninal 
O  of  O3  bonded  to  the  positively  charged  C  atom  of  the  ion, 
followed  by  rapid  02  ejection  and  subsequent  rearrangements 
to  give  the  final  products.  HF  geometry  optimizations  at  the 
6-3 1  G(d)  level  for  POxCl/03"  complexes  show  that  stable 
structures  can  be  found  having  either  the  ozonide-type  or  the 
alkyl  cation— O3  type  structure  for  P02C1",  POCl2",  and 
POCI3".  More  detailed  calculations  of  the  O3  reaction  potential 
energy  surface  and  detailed  reaction  paths  for  the  various 
reactant  ions  are  beyond  the  scope  of  this  work. 

Inspection  of  Table  4  reveals  that  there  is  no  correlation 
between  the  electron  density  on  the  constituent  atoms  of  the 
POxClv”  anions  and  the  reactivity.  This  observation  is  similar 
to  the  computational  results  of  Langeland  and  Werstiuk  regard¬ 
ing  the  reactivity  of  several  neutral  tertiary  phosphite  esters  with 
O3.  In  that  work  using  density  functional  theory  (DFT)  and 
atoms-in-molecules  (AIM)  methods,  no  correlation  was  found 
between  the  size  of  the  activation  barriers  to  produce  phosphates 
and  *02  and  the  atomic  charges  on  the  phosphorus  atom  or 
oxygen  atoms  that  make  up  the  four-membered  ring  of  the 
phosphite  ozonide  intermediate.37 

It  is  interesting  to  note  that  only  the  reaction  of  POCl3"  with 
O3  undergoes  charge  transfer.  Of  the  four  reactants,  only  POCI3" 
has  an  electron  binding  energy  (eBE)  that  is  considerably  less 
than  the  EA  of  O3  (see  Table  1).  The  eBE  of  P02C1"  is  about 
equal  to  that  of  O3,  and  the  charge-transfer  reaction  is  thermo¬ 
neutral  but  is  not  observed.  Not  surprisingly,  the  more  reactive 
ions,  P02C1"  and  POCI3",  have  lower  eBEs  than  the  less  labile 
POCl2".  The  unreactive  P02C12"  ion  has  the  highest  eBE.  A 
similar  observation  was  made  for  the  SOtF;r  species.38  In 


9124-  J.  Phys.  Chem.  A,  Vol.  108,  No.  42,  2004 


Fernandez  et  al. 


addition,  the  reactivity  can  be  explained  by  observing  the  spin 
multiplicities  of  the  anion  reactants  as  shown  in  Table  1  and 
eqs  1  —4.  The  more  reactive  ions  are  doublets,  and  the  less  labile 
species  are  singlets.  Even  though  examples  of  fast  ion— molecule 
reactions  that  do  not  conserve  spin  have  been  observed;39-42  it 
is  apparent  that  both  the  bond  strength  and  spin  multiplicity 
can  affect  reactivity.  Comparable  effects  of  spin  multiplicity 
on  reactivity  were  also  observed  in  an  earlier  study  of  SOJ-)- 
reactions  with  O3,43  where  another  common  feature  of  the 
SOJPjT  and  POxClv-  anions  is  their  lack  of  reactivity  with  CL. 
Finally,  as  noted  above,  there  seems  to  be  no  correlation  between 
reactivity  of  the  anions  and  the  charges  on  the  constituent  atoms. 

Conclusions 

Rate  constants  and  product  branching  fractions  for  the  gas- 
phase  reactions  of  O2  and  O3  with  the  anions  PO2CI-,  POCI3-, 
POCI2-,  and  PO2CI2-  were  measured  from  163  to  400  K  at  0.4 
Torr  in  a  SIFT.  Only  PO2CI-  reacts  with  O2  to  a  measurable 
extent,  while  O3  reacts  with  all  of  the  ions  except  PO2CI2-. 
Stronger  reactivity  is  favored  by  the  following:  (a)  the  weaker 
bond  strength  of  O3  relative  to  O2,  (b)  the  decreasing  eBE  of 
the  more  reactive  anions,  and  (c)  the  higher  spin-multiplicity 
of  the  more  reactive  anions;  i.e.,  doublets  react  faster  than  the 
singlet  species.  There  is  no  correlation  between  reactivity  of 
the  anions  and  electron  density  of  the  individual  atoms.  The 
reactions  of  PCbCl-  with  O2  and  O3  yield  primarily  PO3-  and 
P04“.  The  reaction  of  03  with  POCl2-  yields  mostly  Cl-  and 
PO2CI2-,  while  the  POCI3-  reaction  with  O3  yields  mostly  03- 
and  PO2CI2-.  G3  calculations  were  performed  to  obtain 
optimized  structures,  energies,  and  electron  affinities  for  the 
reactant  and  product  species  and  to  calculate  the  thermochem¬ 
istry  for  the  observed  reaction  channels.  NPA  calculations  show 
that  the  extra  electron  density  added  when  forming  POxCl/" 
anions  is  distributed  primarily  on  the  Cl  atom(s). 

The  results  presented  above  could  have  implications  for  the 
technologies  mentioned  in  the  Introduction.  The  metaphosphate 
anion,  PO3-,  is  interesting  because  it  is  very  stable  and  it  has 
been  detected  in  combustion  environments44  while  the  neutral 
PO3  radical  is  associated  with  phosphorus  impurities  in  coal- 
fired,  magnetohydrodynamic  electric  generators.7’19,44  Impurity 
ions  tend  to  decrease  the  efficiency  of  power  generation  by 
dilution  of  the  plasma  and  by  radiative  loss.  The  present  work 
shows  that  the  reaction  of  PO2CI-  with  O2  and  O3  yield 
significant  amounts  of  PO3-.  However,  the  ozone  is  more  than 
an  order  of  magnitude  more  reactive  over  the  163—400  K  range. 
This  kinetic  information  may  be  of  use  in  the  modeling  of 
impurity  formation  in  plasmas.  PO2CI,  for  example,  has  a 
modest  EA  and  can  be  formed  by  the  reaction  of  neutral  POC1, 
POCI3,  and  PCI3  with  02.45 

The  PO2  moiety  is  thought  to  be  a  catalyst  in  the  recombina¬ 
tion  of  H,  O,  and  OH  flame  species,  which  can  lead  to  an 
enhancement  of  thrust  in  the  supersonic  expansion  in  scramjet 
exhausts.8-12  Apparently,  the  present  reactions  studied  are  not 
capable  of  directly  producing  PO2,  as  no  energetically  allowed 
channels  exist.  However,  reaction  4  has  a  spin-allowed  channel 
that  gives  PO2,  but  it  is  endothermic  by  about  0.5  eV. 

It  is  possible  that  reactions  1  —4  produce  the  neutral  species 
Cl  and  CIO,  which  are  known  to  participate  in  the  catalytic 
destruction  of  stratospheric  ozone.46  This  potentially  negative 
environmental  impact  is  important  when  considering  the  use 
of  oxyphosphorus  additives  for  use  in  combustion  systems. 
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